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Summar 


Current  work  on  this  program  involves  two  projects.  In  the  fir*t  of 
these,  the  effects  of  atmospheric  turbulence  on  target  irradiance  usxng  a 
near-field,  wander-tracking  laser  transmitter  have  been  investigated  in  detail. 
A  unified  analytical  and  phenomenological  treatment  of  the  mean  irradiance 
and  its  fluctuations  is  presented,  with  supporting  experimental  data.  The 
advantages  to  be  gained  through  the  use  of  wander  cancellation  are  seen  to  be 
substantial,  and  in  many  cases  larger  than  predicted  by  the  basic  theory. 

The  discussion  includes  the  observability  of  transmitter-aperture-smoothing 
effects  on  scintillation,  and  the  implications  of  empirical  probability  dis¬ 
tributions  and  power  spectra  for  irradiance-fading  and  wander-tracking  sig¬ 
nals.  Ongoing  analytica  and  multiwavelcngth,  multipathlength  experimental 
efforts  are  expected  to  result  in  a  complete  understanding  of  these  phenomena. 

Regarding  the  second  topic,  a  brief  description  is  also  gi^en  of  recent 
progress  in  the  investigation  of  the  short-term  statistics  of  turbulence  and 
scintillations.  Theoretical  and  experimental  results  have  been  obtained  for 
the  prediction  of  confidence  intervals  or  data  spread  in  such  measurements, 
and  the  effects  of  averaging  times.  In  addition,  a  computer  simulation  tech¬ 
nique  has  been  formulated  for  the  generation  of  an  ensemble  of  instantaneous, 
short  propagation  paths  through  turbulence.  The  results  of  these  efforts  will 
enable  the  prediction  of  limits  on  short-term  scintillation  effects,  and  the 
probabilistic  description  of  such  parameters  as  target  irradiance  given 
apriori  optical/ inf rared  information  on  the  instantaneous  propagation  path. 
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I.  INTRODUCTION 


This  report  describes  progress  on  two  projects.  In  the  first  (Section  II), 
a  detailed  analytical  and  experimental  discussion  is  given  on  target  irradi- 
ance  effects  caused  by  the  action  of  atmospheric  turbulence  on  an  extended, 
wander-tracking  laser  beam.  Ongoing  and  future  activities  are  also  described. 
In  the  second  (Section  III),  a  brief  review  is  given  of  recent  progress  on 
the  problem  of  short-term  propagation  statistics.  The  latter  topic  will  be 
covered  in  detail  in  the  next  technical  report  on  this  program. 


11  *  TURBULENCE  EFFECTS  ON  TARGET  IRRADIANCE  FOR  a  WANDER-TRACKING^  FINITE- 

APERTURE  LASERTRANSMIT'l  ZR  “  ~ 


Fin i  t  e -Aperture  Las er  Transmitter 

Since  the  preceding  report  was  written,  substantial  progress  has  been 
made  in  completing  the  theoretical  problem  of  finite-transmitter  effects, 
including  both  the  phenomenological  and  complete  analytical  descriptions, 
and  their  interrelation.  In  addition,  experimental  data  have  been  obtained 
which  support  these  treatments.  In  this  section,  the  theoretical  treatment 
will  be  described  in  detail,  and  the  experimental  results  given. 

A.  Theoretical  Considerations 

As  described  in  our  previous  reports  on  this  worK,1,2  we  have  taken 
a  dual  theoretical  approach  to  the  problem  of  finite  transmitter  effects  (with 
wander  tracking)  on  target  irradiance.  First,  we  have  approximated  the 
effects  with  a  unified  phenomenological  or  physical  description,  and  second, 
we  have  expanded  and  clarified  the  complete  analytical  treatment  based  on 
the  successful  Huygens-Fresnei  approach3*^  and  reciprocity.3  Although  there 
remain  unsolved  details,  the  overall  topic  is  now  considered  to  be  well 
understood,  as  described  below. 

1.  Mean  Irradiance 

The  general  expression  for  the  mean  irradiance  is  given  by3,/*’6 

ik - 

-  -  /  k  \2  2-  -  P'P 
1  (P)  "l  1™  )  /d  P  M(P)  e 


X  /  d2R 


4+iHr-iK 


ik  —  — 

~  p.  R 


(1) 


1. 

J.  R. 

pher ic 

1. 

J.  R. 

pher ic 

3. 

R.  F.  ] 

4. 

H.  T.  ’ 

5. 

D.  L.  ] 

6. 

H.  T.  ' 

°Pt.  10,  2771  (1971) 


Opt.  10,  1652  (1971), 


-4- 


where 


M 

and  U 


transverse  position  vector  for  target  point  under  consideration 
optical  wavenumber 
path  length 

sum  ard  difference  coordinates  in  the  transmitter  aperture, 
respectively 

atmospheric  modulation  transfer  function 
complex  transmitter  amplitude  distribution. 


We  immediately  specialize  the  transmitter  to  the  case  of  a  gaussian  beam  with 
nominal  radius  (a)  and  output  focal  length  or  radius-of-curvature  f: 


U(r) 


U  e 
o 


1 

2 


2 

r 


(2) 


The  second  integral  in  Eq.  vl)  then  becomes 


(3) 


We  note  immediately  that  the  (R>p)  exponential  term  vanishes  in  the  case 
of  a  focused  beam  (f=z),  including  focus  on  infinity  (collimation) .  Thus  the 
expression  may  be  immediately  compared  with  that  for  the  reciprocal  case  of  a 
heterodyne  receiver,  in  which  a  plane  wave  (or  infinite  focus)  has  generally 
been  assumed. ^ 

We  now  use  Eq.  (3)  with  (1)  to  write 


7.  D.  L.  Fried,  Proc.  IEEE.  55,  57-67,  January  1967. 
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We  let  d  p  =  pdp  d<f>  and  integrate  over  <J>  to  obtain  a  single  integral  equation: 


I  (P)  = 


kV 

2z2 


CO 

|Uol 2  5  dp  p  M(p)  J0  (z  PP) 


2 

_  P_ 
4 


a  z 


(5) 


This  consUtutes  a  generalization  of  Eq.  (A3)  in  Ref.  6. 

We  now  consider  the  special  case  of  interest  here,  in  which  the  target 
point  is  on-axis  (p=0)  and  the  transmitter  is  focused  on  the  target  (f=z): 

*  *  |U0|2  *  dp  p  M(p>  e"P  •  (6) 

2z  o 

It  may  be  pointed  out  that  the  dependence  on  the  Fresnel  number  (ka2/^)1^2  = 

8  has  disappeared  from  the  integral:  the  pertinent  independent  variable  will 
simply  be  (a/pQ),  where  pQ  is  the  coherence  scale  owing  to  turbulence.  This 
means  that  basic  conclusions  drawn  from  earlier  analyses7  of  the  reciprocal, 
heterodyne  receiver  case  will  be  valid  for  a  focused  or  collimated  trans¬ 
mitter  with  f-z,  over  an  arbitrary  (nonuniform)  path;  it  is  only  necessary 
to  know  M(p)  for  the  path.  If  we  define  the  transmitter  diameter  D=2a,  and 
also  define  xHp/D,  we  can  write  Eq.  (6)  in  the  notation  of  Ref.  6: 
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I  *  2B^  |UqJ  /  dx  x  M(Dx)  e  x  .  (7) 

o 

For  the  collimated  case  with  z<°°,  the  exponent  in  the  integrand  is  multiplied 
2 

by  (1  +6  ).  Finally,  for  <the  vacuum  or  turbulence-free  case,  we  let  M*1  and 
find 


I 
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oo  2 
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If  wander-tracking  is  not  employed,  we  have  the  so-called  "long-term 
case",  for  which  the  atmospheric  modulation  transfer  function  is3’^’6 


M(p) 


M(Dx) 


(9) 


where  is  the  turbulence-induced  coherence  scale.  The  complete  expression 
for  the  irradiance  is  thus 


/ 


o 


(Eq.  A8  of  Ref.  6). 


(10) 
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For  the  wander-tracking  case,  the  "short-term  modulation  transfer  func¬ 
tion  (M^)"  is  required,  which  is  dependent  on  the  optical  configuration,  and 
which  is  related  to  the  reciprocal  problem  of  wavefront-tilt-tracking  in  a 

8  9 

heterodyne  receiver  or  short-exposure  resolution  in  an  imaging  system. 

Unfortunately,  this  problem  has  not  been  adequately  solved.  In  Ref.  6,  a 

8  9 

result  for  M  is  used  which  was  derived  by  Fried  *  for  a  uniformly-weighted, 
truncated  aperture: 


If  this  is  substituted  in  Eq.  (7),  it  is  seen  that,  for  D  much  greater  than 
P  t  the  integral  diverges;  i.e.,  the  equivalent  structure  function  in  the 
exponent  becomes  negative. 

The  basic  problem  here  is  that,  in  the  derivation  of  Eq.  (11),  a  certain 
assumption  was  made  regarding  the  independence  of  the  wavefront  tilt  and 
higher-order  phase  distortions  over  the  aperture.  If  this  assumption  is 
dropped,  higher-order  terms  appear  in  the  exponent  of  (11),  which  prevent  the 
exponent  from  becoming  positive  and  hence  diverging  the  integral.10  Although 
the  correction  may  be  small  for  the  case  of  a  truncated,  uniform  aperture,  it 
is  mathematically  and  physically  crucial  for  the  inf inite-gaussian  beam.  A 
complete  expression  for  the  latter  aperture-distribution  is  under  investiga¬ 
tion  by  R.  Lutomirski. 11 

Unfortunately,  the  results  of  Ref.  6  for  the  wander-tracked  problem  must 
be  entirely  discounted.  We  summarize  the  reasons  as  follows: 

(1)  The  concept  of  a  "si  ort-term  coherence  scale"  is  not  justified, 
since  the  1/e  radius  of  M,,^  is  not  meaningfully  related  to  a 
"short-term  beam  spread".  1 


8.  D.  L.  Fried,  "Effects  of  Atmospheric  Turbulence  on  Static  and  Tracking 
Optical  Heterodyne  Receivers/Average  Antenna  Gain  and  Antenna  Gain  Varia¬ 
tion",  Technical  Report  No.  TR-027,  Optical  Science  Consultants,  Aug. 1971. 

9.  D.  L.  Fried,  J.  Opt.  Soc.  Am.  56,  1372  (1966). 

10.  R.  Lutomirski,  "The  Tilt-Corrected  MTF",  Preliminary  Report,  Pacific- 
Sierra  Research  Corp.,  Santa  Monica,  California,  March  1974. 

11.  D.  L.  Fried  and  R.  Lutomirski,  private  communications. 
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(2)  The  function  MgT  used  was  inappropriate  for  the  gaussian  case 
considered . 

(3)  A  simple  MacClaurin  expansion  will  show  that  Eq.  (12)  of  Ref. 

6  represents  a  good  approximation  to  Eq.  (11)  only  over  a 
very  limited  and  uninteresting  range  of  Pq/D,  even  though  the 
figures  in  the  paper  imply  a  much  wider  range  of  validity. 
Specifically,  Pq/D  must  be  smaller  than  0.03  if  the  second- 
order  term  is  to  be  less  than  10%  of  the  first-order  term 
used  in  Eq.  (12)  of  that  paper. 

Although  we  do  not  yet  have  a  correct  M^,  we  point  out  that  the  gaus- 
sian-beam  case  will  not  be  significantly  different  from  that  for  a  uniformly- 
weighted  aperture,  and  that  Fried's  corresponding  result  is  expected  to  be  a 
good  approximation.  This  is  further  supported  in  the  discussions  below  of 
the  phenomenological  model  and  the  experimental  data. 

Wc  now  distinguish  two  cases.  Although  Eq.  (10)  may  be  plotted  uni¬ 
versally  as  a  function  of  (D/pq),  the  normalization  by  Tq,  the  turbulence- 
free  or  D/po+0  value,  implies  that  B2  -  kV/l6z2  remains  constant  (note 


that  Iq  -  lujV  ^ 
16z 


Transmitter  Power 


1  2n2 

R  D 

X  Hence,  the  implied 


indepe  ent  variable  is  Pq,  i.e.,  the  turbulence  strength.  As  an  alternative 
normalizer,  let  us  consider  the  limit  in  the  case  D/p  -*•  «.  We  define 
o  -  D/pq,  and  consider  the  limit  of  the  integral  in  Eq.  (10): 


"  2  5/3 

lim  / dx  x  e  X  e  ^ax^ 
a-*®  o 
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Hence  we  may  write 
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D/po 
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^  Transmitter  Power  X 


(14) 


In  particular,  this  justifies  the  viewpoint  of  a  diffraction  spread  from  a 
scale  length  (pQ)»  i.e.,  proportional  to  (1/kp^)2,  and  it  furthermore  pro¬ 
vides  a  definite  numerical  constant  for  the  asymptote  in  Fig.  la  of  Ref.  1. 
This  will  be  further  discussed  below.  The  expression  corresponding  to  Eq. 
(10)  is 
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(15) 


and  the  implied  independent  variable  is  D,  such  that  the  normalizer  I  is 

00 

constant  with  constant  transmitter  power. 

It  is  instructive  to  compare  the  results  (8)  and  (14)  with  corresponding 
results  for  the  reciprocal  cases  in  the  literature— noting  that  the  latter  are 
for  a  uniform,  truncated  aperture.  Anticipating  a  simple  result,  let  us  iden¬ 
tify  the  effective  diameter  D'  of  an  equivalent  uniform  truncated  aperture 
with  our  gaussian  D  as  follows: 


D'  -  2D  -  4a 


(16) 


We  also  note  that  parameter  r  as  utilized  in  Refs.  7-9  is  related  to  p  bv 

o  o  J 


r  *  2.0986  p 
o  o 


(17) 


We  then  state  the  results  for  three  cases: 
(1)  Gaussian  transmitter  (Eqs.  8,  14): 
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1.1 


(18a) 
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(2)  Uniform,  truncated  image  sy-stem: 
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(3)  Uniform,  truncated  heterodyne  receiver: 
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The  uniformity  of  the  large-to-small  aperture  ratios  in  all  three 

constitutes  the  basis  for  Fq.  (16),  and  we  note  that  i.i  p  ^ 
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Due  to  compatible  dimensionality,  Eqs.  (18a)  and  (18b)  provide  a  meaningful 
numerical  comparison  of  the  gaussian  and  uniform,  truncated  cases  respectively. 

The  numberical  evaluation  of  Eqs.  (10)  and  (15)  is  shown  in  Fig.  1.  Also 
shown  is  the  corresponding  result  for  the  uniform-aperture  imaging  system;^ 

Eqs.  (16,17)  have  again  been  utilized,  the  former  representing  a  shift  to  the 
left  by  one  octave. 

An  analysis  of  the  long-term  and  wander-tracked  cases  has  appeared  in  the 
Russian  literature. 12  In  this  derivation,  Kon  treats  a  focused,  gaussian 
beam,  and  utilizes  an  expression  for  mean  irradiance  which  is  essentially 
identical  to  that  of  the  Huygens-Fresnel  formulation.  For  his  wander-tracked 
results,  he  defines  the  beam  centroid  in  a  straightforward  manner  and  calcu¬ 
lates  I  following  this  centroid.  His  results  are  given  in  terms  of  1/1  vs. 

the  phase  structure  function  (2a)  in  degrees;  the  relationship  to  our 
abscissa  is 


/VD) 


(degrees) 
307.1  ' 


3/5 


(19) 


Unfortunately,  his  long-term  results  disagree  substantially  with  those  of 
Eq.  (15)  and  Ref.  9,  and  his  tracking  results  are  similarly  in  doubt.  In  par 
ticular ,  his  results  show  ^90%  fall-off  in  ILT  for  -  0.25  (cf.  Fig.  1). 

As  will  be  seen  below,  experimental  wander-tracked  results  to  date  agree 
sufficiently  with  the  predictions  of  Ref.  8  to  indicate  that  the  gaussian  and 
uniform-truncated  cases  are  not  substantially  different. 

We  return  now  to  the  phenomenological  model  discussed  in  Refs.  1,2.  If 
we  represent  the  total  solid-angular  bearaspread  as  consisting  of  terms  owing 
to  diffraction  [„  1/k2,2]  .  wander ^  (2a)/k2a2  ^1/3  ^5/3^ -1!  ,  and 

beam  breakup/  1/R2p  2\  ,  we  can  write  the  mean  irradiance  as 


k^» 


1. 

3  +  C. 


.2  2 

k  po 


12.  A.  I.  Kon, 


Izv.  VUZ  Radiofizika  n,  61  (1970). 
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(20) 


14 


Figure  1.  Theoretical 
Ref.  9. 


where  the  coefficients  remain  to  be  determined.  We  are  implicitly  defin¬ 

ing  our  angular  sizes  in  terms  of  their  effect  on  T,  which  removes  the  compli¬ 
cation  of  alternative  angular-  beamwidth  definitions.  We  note  that  for  our 
gaussian  beam  (see  Eq.  2), 


I  (r) 


Transmitter 


!Uol 


2.  2 
■r  /a 


where  |U  1 
o 


2 


Transmitted  Power 
2 

ira 


(21) 


Diffraction-Limited  Case  (a to  ■*  0) 

.  - -  c 

From  basic  optics,  we  have  at  the  target  (focal  point): 


I 


,2  * 

k  a 


2 

z 


z2 

where  w  «  =  new  gaussian-beam  radius 

k  a 


(22) 


from  which  we  conclude  that 


=  1 


(23) 


for  a  gaussian  beam. 

Beam-Breakup  Case  (a/p^  00 ) 

From  Eq.  (14),  letting  a  ■  D/2,  we  immediately  have 


C3  =  3.63 


(24) 
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(As  an  aside,  we  note  that,  since  the  transverse  irradiance  distribution  is 
probably  gaussian  for  this  case,  the  gaussian-beam  radius  corresponding  to 


Eq.  (24)  is 


w'  =  ^3.63z2/k2pQ2 


1/2  .) 


We  now  note  that  the  asymptotic  intersection  of  the  two  above  cases 
occurs  at 


a  =  0.525p  ■  r  /4 

o  o 


D  -  1.05p 

o 


(25) 


For  the  uniform,  truncated  aperture  case,  this  breakpoint  occurs  at  the 
fundamental  point  D  *■  r^,  which  again  justifies  our  taking  2D  as  the  equi¬ 
valent  uniform  aperture  diameter  (Eq ,  18). 

Wander  Case 

The  determination  of  the  constant  C ^  is  less  straightforward. 

It  can  be  derived  crom  approximate  physical  arguments,  as  discussed  below, 
or  from  fitting  Eq.  (20)  to  the  numerical  results  from  Eq.  (10).  To  facili¬ 
tate  this  comparison,  we  rewrite  (20)  as 


(26a) 


which  correspond  to  Eqs.  (10)  and  (15)  respectively.  The  best  fit  occurs  if 
C 2  is  obtained  from  the  breakpoint  described  above,  and  results  in  the  value 
(C2  =  0.423).  The  resulting  curves  are  essentially  indistinguishable  from 
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those  of  Fig.  1.  An  approximate  physiccl  analysis  from  Titterton,^  ^sing 
Kon's  mean-square  wander  angle  (0. 76Cn2za~1/,3) ,  results  in  C2  *  1.39  and 
a  less  satisfactory  fit  to  the  curve. 

It  should  he  pointed  out,  however,  that  the  values  obtained  from  Eqs. 

(26a, b)  are  insensitive  to  C 2<  In  fact,  a  good  approximation  is  to  let  C2  * 

0  (see  Fig.  2);  this  suggests  that  most  of  the  wander  effect  is  contained  in 

2  2 

the  1/k  term  of  Eq.  (20). 

Since  we  wish  to  use  this  phenomenological  description  for  the  wander- 
tracked  case,  and  most  of  the  wander  is  contained  in  this  p  “2  term,  we  postu¬ 
late  that  the  proper  coefficient  for  the  tracked  case  is  negative .  Al¬ 
though  we  do  not  yet  have  a  reliable  H aygens-Fresnel  result  for  the  gaassian- 
beam  case,  as  indicated  above  we  believe  that  the  truncated,  uniform  aperture 
results  will  be  substantially  similar.  In  Fig.  3,  we  plot  Eq.  (26a, b)  with 
C2  -  -3.44,  and  the  similar  curves  deduced  from  Ref.  8.  This  value  was  chosen 

so  that  the  maximum  from  Eq.  (26b)  would  agree  with  that  from  Ref.  8,  which 

predicts  a  gain  of  3.4  (5.3  dB)  in  target  irradiance  or  10.6  dB  in  receiver 
photocurrent  over  that  for  a  very  large  transmitter.  The  peak  of  the  latter 
curve  occurs  at  D/pq  '  4  or  D/rQ  *  1-9;  this  is  consistent  with  the  antenna 
gain  realized,  and  indicates  that  half  the  "effective"  aperture  of  2D  is  used 
coherently  at  this  peak.  As  will  be  seen  below,  the  tracking  improvetm  it  in 
irradiance  fading  is  also  optimized  for  3  'v  D/p  'v  5.  We  point  out  that  the 
curves  of  Fig.  3  are  generally  consistent  with  those  of  Figs,  la,  lb  in  Ref.l. 

Finally,  we  plot  the  ratio  of  the  mean  irradiance  with  tracking  (1^,)  to 
that  without  tracking  (1^,)  in  Fig.  4,  again  using  the  u.iiform-aperture 
results  of  Ref.  8.  In  lieu  of  a  result  for  the  gaussian  beam  case,  this 
curve  will  be  used  for  comparison  with  experimental  results  below.  The  maxi¬ 
mum  value  is  predicted  to  be  6.2  dB.  It  may  be  pointed  out  that  the  peak  in 
this  curve  suggests  a  possibility  for  instrumenting  a  "po-meter",  utilizing  a 
point  source  and  an  image-tracking,  variable  aperture  test  receiver  at  the 
target  and  transmitter  respectively. 

*In  the  remainder  of  this  report,  we  identify  the  wander-tracked  mean-irradi- 
ance  and  variance  with  the  subscript  T,  and  the  corresponding  non-tracking 
values  with  the  subscript  NT,  where  the  latter  is  identical  to  the  "long¬ 
term  case.  Also,  in  plotting  the  results  from  Ref.  8,  we  have  interpreted 
the  diameter  in  that  reference  as  being  twice  our  D,  as  explained  above,  and 

have  corrected  for  the  use  of  p  in  the  place  of  r  . 

_ _ _  O  o 

13.  P.  J.  Titterton,  Appl.  Opt.  12,  423  (1973). 
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Figure  3.  Theoretical  mean  irradiance  with  wander-tracking 


2.  Fluctuations  in  Irradiance  (Fading) 

The  problem  of  fading  due  to  beam  wander  and  scintillation  is 
considerably  more  difficult.  In  the  Huygens-Fresnel  formulation,  this  second 
moment  of  irradiance  (fourth  moment  of  amplitude)  involves  an  eight-fold 
integral.  The  basic  expression  is  given  as  Eq.  (42)  in  Ref.  14,  or  equi¬ 
valently,  as  Eq.  (2.15)  in  Ref.  15.  Until  recently,  the  only  attempt  to 

evaluate  this  expression  was  in  Fried's  analysis  of  "atmospheric  modulation 

8  15 

noise"  (coherent  fading)  in  a  heterodyne  receiver;  *  however,  the  approxi¬ 
mations  used  in  these  treatments  lead  to  erroneous  results. 

Interesting  work  on  this  problem  was  recently  reported  \y  Russian 
investigators. I**  Using  the  assumptions  that  the  atmospheric  perti  rbation  of 
the  propagation  Green's  function  is  essentially  a  phase  effect,  and  that  the 
phase  dis_.ibution  is  gaussian,  they  reduced  the  problem  to  a  six-fold  inte¬ 
gral  which  they  then  evaluated  using  Monte  Carlo  techniques.  The  result  will 
be  seen  below  to  agiee  in  many  respects  with  our  phenomenological  description. 

In  particular,  it  predicts  a  v^lue  of  unity  for  the  normalized  irradiance 

2  1 
variance  fa^  )  at  large  values  of  D/ Pq,  which  agrees  with  our  predictions. 

This  indicates  that  the  Huygens-Fresnel  formulation  is  valid  under  conditions 

of  large  phase-distortions. 

The  theoretical  results  of  Ref.  16  can  be  converted  to  our  param¬ 
eters  by  squaring  the  ordinate  and  relating  the  abscissa  as  follows: 


D  /  (Abscissa,  Ref.  16)  \ 
°o\  JT  ) 


The  results  are  shown  in  Fig.  5.  A  major  point  is  that,  for  a  focused  beam, 
the  treatment  of  Ref.  16  results  in  a  curve  which  is  a  universal  function  of 
the  spherical-wave  phase  structure  function  or  D/pq.  We  believe  that  this 
is  not  entirely  correct,  and  that  the  transmitter  size  or  Fresnel  number  is 
explicitly  involved,  as  explained  in  the  phenomenological  description  below. 
This  means  that  different  curves  sr.culd  be  derived  depending  upon  which  param¬ 
eter  (D  or  pQ)  is  varied;  the  results  of  Ref.  16  are  apparently  best  related 
to  the  fixed-D,  variable- p  case. 


14.  H.  T.  Yura,  Applied  Optics  H.,  1399-1406  June  1972. 

15.  D.  L.  Fried,  J.  Quantum  Electr.  QE-3,  213  (1967). 

16.  V.  A.  Banakh,  et  al,  J.  Opt.  Soc.  Am.  64,  516-518,  April  1974. 


zed  variance  of  irradiance  without  wander-tracking,  from  Ref.  16 


The  fact  that  the  assumptions  used  in  Ref.  16  lead  to  inexact  results  is 
clearly  seen  through  inspection  of  Eq.  (17)  of  that  reference  for  the  mean 
irradiance:  it  is  identical  to  the  exact  expression  (Eq.  5  of  the  present 
report,  with  p  -  0)  only  for  the  focused  case.  The  inexactness  is  probably 
related  to  the  dropping  of  the  amplitude  term  in  the  perturbation  function. 

We  intend  to  pursue  the  implications  of  this  interesting  approach  fur¬ 
ther,  and  to  attempt  to  obtain  a  numerical  solution  for  the  full  eight-fold 
expression.  It  will  also  be  of  great  interest  to  incorporate  the  wander¬ 
tracking  results  when  they  are  available.  In  the  latter  case,  for  a  focused 
beam,  the  primary  point  of  interest  is  whether  or  not  the  "transmitter  smooth¬ 
ing  of  scintillations"  is  predicted  by  the  Huygens-Fresnel  formulation,  in 

agreement  with  first-order  theory. ^  ^  This  point  will  be  discussed  further 
below. 

We  now  review  the  phenomenological  description  of  the  fading,  with  some 
added  considerations  over  those  reported  in  Refs.  1,2.  We  discuss  three  inde¬ 
pendent  fading  mechanisms:  wander-fading,  first-order  scintillation,  and 
coherent  fading  (beam-breakup  scintillation). 

Wander  Fading 

Using  Titterton's  development  in  Ref.  13,  we  write  the  normalized 
irradiance  variance  due  to  wander  as 


Wander  4y  +  1 


(28) 


where  y  is  defined  as  the  ratio  of  the  mean-square  wander  angle  to  the 
short-term  beamspread  angle.  Using  Eq.  (20),  we  write  this  as 


17.  A.  Ishimaru,  Radio  Science  295  (1969). 

18.  D.  h.  Fried  and  R.  A.  Schmeltzer,  Applied  Optics  6,  1729-1737,  Oct.  1967. 

19.  J.  R.  Kerr  and  R.  Eiss,  J.  Opt.  Soc.  Am.  62,  682  (1972). 
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(29) 


where  we  interpret  O,  as  the  negative  number  relating  to  the  wander-tracking 
results  discussed  above;  this  coefficient  is  taken  as  the  best  measure  of  the 
effective  wander.  We  then  combine  Eqs.  (28,  29)  to  write 


(30) 


For  the  particular  choice  of  C.,  discussed  above  (-3.44),  this  becomes 


(30a) 


We  note  that  the  wander-fading  term  is  a  universal  function  of  D/p  . 

The  general  behavior  of  Eq.  (30)  can  be  seen  by  approximating  the 
function  in  terms  of  four  log-linear  realms  (cf.  Fig.  21,  Ref.  2).  This  is 
shown  in  Fig.  6.  The  breakpoints  occur  at  values  of  D/p  on  the  order  of 

1/2,  1,  and  100,  with  a  maximum  value  of  Oj2  on  the°order  of  (in  excess 

of)  unity.  wander 
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ig  due  to  beai 


First-Order  Scintilla t ion 


By  first-order  scintillation,  we  mean  that  which  would  be  pre¬ 
dicted  for  a  focused  transmitter  by  the  Rytov  approach  (Refs.  17-19).  This 
includes  the  "transmitter  smoothing  of  scintillations".  It  is  now  clear1’2, 
that  in  order  for  this  effect  to  exist,  it  is  necessary  that  p  >  D>>(z/k)1^2; 
the  requirement  that  Pq  be  greater  than  (z/k)1^2  is  dimensionally  equivalent 
(Eq.  32  below)  to  requiring  that  the  scintillations  from  a  point  source  be 
weak  (unsaturated).  As  explained  below  in  conjunction  with  the  experimental 
data,  it  also  appears  that  the  inequality  requirement  (p  /D>1)  may  in  fact  be 
(Pq/D>>1):  a  small  amount  of  atmospheric  phase  distortion  over  the  aperture 
may  destroy  the  smoothing  effect. 

The  customary  way  to  write  this  fading  term  is1,8,15 

2  b(t>'o  2 

°I  E  e  X  -  1  (31) 

ScintilJation 


where  is  the  log  amplitude  variance  for  a  point  source,  and  B'is  the 
smoothing  factor.  In  writing  this  it  is  assumed  that  the  fluctuations  are 
log  normal.  3'  is  a  function  of  D/(z/k)1/2  or  the  Fresnel  number,  and  for 
large  Fresnel  numbers19  is  proportional  to  D~7/3.  Noting  that 


2 

a 

X 


0.228 


(32) 


and  assuming  that  8'o  2«1,  we  may  write  (31) 

X 


2 

Scintillation 


£  0.912  S' 


-  0.912 


5/3 


(33) 


This  is  obviously  not  a  universal  function  of  D/p  . 

20.  K.  S.  Gochelashvily,  "Focused  Irradiance  Fluctuations  in  a  Turbulent 
Medium",  Optica  Acta  20,  193  (1973). 
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20 


Coherent  Fading 

Coherent  fading  or  beam-breakup  scintillation  corresponds  in  a 

O  1C 

reciprocal  heterodyne  receiver  to  "atmospheric  modulation  noise".  *  It 

becomes  predominant  when  D  >  some  fraction  of  p  .  The  correct  form  for  this 

—  o 

fading  term  is  not  known,  but  from  the  results  of  Ref.  16  and  physical  reason¬ 
ing  we  expect  it  to  be  an  increasing  and  perhaps  nearly  universal  function  of 
D/pq.  However,  contrary  to  Refs.  (8,  15),  we  do  not  expect  it  to  increase 

without  limit.  For  large  values  of  D/p  ,  we  invoke  an  analysis  of  Brown's21 
2  2  ° 

involving  N  ■  D  /pq  oscillators  having  identical  frequencies  but  independent 
mutual  phases;  the  result  is 


o 


2 

Coherent  Fading 


'V/ 


»1 


(34) 


which  approaches  unity  for  large  (D/p  ).  The  expression  is  again  not  a  uni- 

o  j 

versal  function  of  (D/p  ) ,  since  it  involves  a  as  given  in  Eq.  (32).  For 
2  °  X 

o  <<  1  it  is  given  by 
X 


'V 

'V 


1  - 


fe 


Coherent  Fading 


»1,  a  «1 


) 


1  +  4o 


»o2 


(34a) 


2 

For  values  of  a  exceeding  0.17,  the  expression  (34)  Indicates  that 
2  x 

°T  .  ,  „  exceeds  unity  for  sufficiently  large  D/p  .  The  general 

Coherent  fading  o 

behavior  of  the  coherent  fading  term  is  shown  in  Fig.  7,  where  the  hump  is 

2 

seen  to  be  significant  as  o  approaches  the  saturation  value  of  0.6. 

X 


21.  W.  P.  Brown,  Jr.,  "Research  in  Interaction  of  Coherent  Light  with  Solids 
and  with  Turbulent  Atmospheres",  Research  Report,  Hughes  Research 
Laboratories,  Malibu,  Calif.,  May  1972. 
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Normalized  irradiance  variance  due  to  coherent  fading  for  curve 


Combined  Fading 

In  order  to  consider  combined  fading,  we  add  Eqs.  (30),  (31  or  33), 
and  the  curve  of  Fig.  7.  Since  the  latter  two  components  are  not  universal 
functions  of  D/p  ,  we  must  distinguish  whether  the  aperture  size  or  the 

2 

coherence  scale  are  being  varied.  In  the  former  case,  we  must  specify  o 

1/2  X 

or  pQ/(z/k)  ;  we  then  obtain  curves  such  as  that  shown  in  Fig.  8.  In  the 

latter  case,  we  must  specify  the  Fresnel  number,  and  we  obtain  curves  such  as 

that  shown  in  Fig.  9;  this  result  with  wander  effects  is  in  general  agreement 

with  that  of  Fig.  5. 

With  the  exception  of  the  exact  behavior  of  coherent  fading,  and 

the  knowledge  of  the  best  value  for  C_,  sufficient  information  is  given  here 

L  2  2 

to  construct  a  physically-based  prediction  of  vs.D  (given  o  )  or  vs. 

pQ  (given  Fresnel  number).  Furthermore,  the  predicted  behavior  with 

wander  tracking  is  readily  determined  by  dropping  the  contribution  from  that 

mechanism.  In  particular,  the  curve  of  Fig.  9  will  be  shown  later  to  agree 

well  with  both  the  tracking  and  non-tracking  experimental  data  on  fading. 

Also,  the  curve  of  Fig.  8  suggests  an  alternative  "p  -meter"  to  that  mentioned 

o 

above;  a  non- image-tracking  receiver  would  be  utilized  to  determine  the  aper¬ 
ture  size  D  for  maximum  fluctuations  in  the  photocurrent. 

Finally,  we  emphasize  that,  in  the  presence  of  wander  tracking, 
there  is  definitely  an  optimum  aperture  size  from  the  standpoint  of  fading  as 
well  as  mean  irradiance,  and  too  large  an  aperture  can  have  significantly 
deleterious  effects  on  fading  performance. 

3.  Other  Aspects 

Generalized  Focus  Conditions;  Nonuniform  Turbulence  (Vertical 
Propagation) 

For  the  case  of  a  general  turbulence  profile  over  the  path,  the 
atmospheric  effect  on  the  mean  irradiance  is  still  represented  by  Eq.  (9), 
where 


1.45  k2 


-  3/5 


(35) 


and  the  integration  is  from  the  target  to  the  transmitter.  The  degrading 
effect  of  a  given  level  of  turbulence  on  transmitter  coherence ^ p  ^  is 

most  heavily  weighted  towards  the  transmitter  end.  Hence,  an  uplink  beam 
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irradiance  vs.  transmitter  diameter,  with  and  without  wander-cancel lat ion 


Normalized  variance  of  irradiance  vs.  coherence  scale,  with  and  without  wander  cancellation 


tends  to  be  incoherent  or  broken  up,  while  a  downlink  beam  remains  intact. 

In  any  event,  the  nonunifirm  distr  bution  of  turbulence  over  the  path  is 

entirely  accounted  for  in  the  parameter  p  . 

o 

Let  us  now  consider  the  effect  of  an  arbitrary  beam  focus.  We  write 
the  on-axis  mean  irradiance  from  Eq.  (!>)  as 


LT 


00  -(p/p  )  5/3  -  p2/4 

\+lc2.2  (i-i)2 

/  dp  p  e  V  07  e 

2  V  z  f  / 

a 

/  dp  p  e 


-(p/%) 


5/3  2 


-  P  ia’2 


(36) 


where 


a  = 


V(M)  - 


2  2 
a 


1 

2 


(37) 


is  the  generalized  beam  parameter.  The  defocusing  term  will  reduce  and 
wholly  determine  a'  unless  a  precise  focus  adjustment  (z  *  f)  is  used.  In 
any  event,  if  pQ  >>  a',  the  beam  term  predominates  in  Eq.  (36);  this  will 
appiy  to  a  downlink  in  particular.  If  p^  <<  a',  the  atmospheric  term  pre¬ 
dominates  regardless  of  (reasonable)  beam  adjustment;  this  applies  to  a 

large-aperture  uplink.  In  any  far-field  case,  the  actual  focus  adjustment 

2  2  2  o 

is  not  important  as  long  as  k  a  /f  is  smaller  than  1/a  (Eq.  37). 

An  important  point  is  that,  once  a'  is  determined  from  Eq.  (37)  and 
PQ  known,  the  mean  irradiance  is  determined  from  the  same  curve  as  in  the 
horizontal,  focused  case.  Hence  this  formulation  greatly  simplifies  the 
treatment  of  the  nonuniform,  general-focus  case. 

For  the  wander-tracked  case,  there  will  appear  an  optical-system- 
dependent  correction  to  the  atmospheric  term  in  Eq.  (36).  This  will  be 
significant  only  for  a  on  the  order  of  p^,  which  can  be  achieved  with  an 
uplink  beam  having  a  limited  aperture  size. 
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Similar  reasonings  can  be  applied  to  the  fading.  It  appears  that  the 
results  will  again  indicate  that  the  turbulence  profile  is  entirely  contained 
in  pQ.  However,  for  a  general,  defocused  beam,  the  parametric  dependencies 
are  complicated.  From  the  first-order  theory,22  which  implies  wander-track¬ 
ing,  we  expect  significant  transmitter  smoothing  of  scintillations  if  the 
aperture  is  smaller  than  pq  and  the  turbulent  region  is  confined  to  the  near 
field  of  tne  transmitter.  This  has  not  yet  been  demonstrated  from  the  Huy- 
gens-Fresnel  equations. 

Outer  Scale  Effects 

We  have  so  far  in  this  discussion  neglected  effects  relating  to 
the  limitation  on  the  extent  of  the  inertial  subrange  of  turbulence.  Luto- 
mirski  and  Yura  have  considered  the  effects  of  the  outer  scale, 2^,2^and  in 
particular  they  derived  first-crder  corrections  to  such  quantities  as  turbu¬ 
lence-induced  beamspread  and  coherence  scale  pq,  where  "corrections"  refer 

to  the  use  of  the  modified  von  Karmann  in  place  of  the  (nontruncated)  Kol- 
23 

mogorov  spectrum. 

The  critical  parameter  in  their  development  is  the  ratio  (z/z  ), 

c 

where  z c  is  a  critical  pathlength  dependent  on  the  outer  scale  (L  ),  and  the 
ratio  is  given  by 


In  this  expression,  pQ  is  the  uncorrected  or  Kolmogorov  value.  Hence,  the 

outer  scale  effects  are  implied  to  be  dependent  only  on  L  and  p  ,  and  not 

explicitly  related  to  z,  C  2,  k,  or  D. 

n 

The  percentage  correction  of  pQ  is  shown  as  approximately  182, 
252,  and  502  for  z/z ^  *  100,  10,  and  2  respectively.  The  corrections  are 
such  as  to  increase  pQ,  since  the  von  Karmann  spectrum  has  less  energy  at 
small  wavenumbers  than  does  the  Kolmogorov. 

22.  J.  R.  Kerr  and  J.  R.  Dunphy,  J.  Opt.  Soc.  Am.  63,  1  (1973). 

23.  R.  F.  Lutomirski  and  H.  T.  Yura,  J.  Opt.  Soc.  Am.  61,  482-487, 

April  1971.  — 

24.  H.  T.  Yura,  J.  Opt.  Soc.  Am.  83,  107-109,  January  1973. 
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To  avoid  outer  scale  effects  in  experiments,  it  is  apparently  necessary 
to  maximize  z/z £.  For  a  given  outer  scale,  this  indicates  small  values  of 
PQ.  However,  it  may  be  desirable  to  work  in  the  realm  Pq>  D  while  maintain¬ 
ing  a  large  Fresnel  number^  -  D  (k/z)1/2J  .  We  therefore  write  Eq.  (38)  as 


z_ 

z 

c 


15.3  L 


5/6 


(39) 


which  as  expected  shows  that  small  wavelengths  and  pathlengths  reduce  outer- 
scale  effects.  Conversely,  in  order  to  perform  experiments  on  outer  scale 
effects,  large  wavelengths  and/or  pathlengths  are  required. 

For  the  1.6  km  pathlength  and  6328  X  wavelength  utilized  in  the  experi¬ 
ments  to  be  described  below,  and  with  Lq  ■=  lm,  we  have  z/zq  >  100  in  all 
cases.  At  10.6  p,  with  pq/D  =  5,  6  -  5,  and  Lq  *  lm,  we  have  z/z^  =  100  at 
a  pathlength  of  z  =  100  m.  In  any  event,  it  should  be  noted  that^his  treat¬ 
ment  predicts  a  surprisingly  large  outer  scale  effect,  even  for  z  »  z  . 

25  c 

In  recent  work.  Greenwood  and  co-workers  have  empirically  determined 

a  small-wavenumber  turbulence  spectrum  which  is  intermediate  between  the 
Kolmogorov  and  von  Karmann.*  It  would  be  useful  to  recompute  structure  func¬ 
tions  and  pQ  based  on  this  spectrum.  In  particular,  the  expected  effect  on 
beam  wander  should  be  computed. 

As  discussed  below,  our  data  (and  those  of  others)  clearly  indicate  a 
larger  component  of  wander  than  predicted  by  our  expressions,  especially  at 
large  values  of  D/pQ,  and  this  is  thought  to  be  related  to  large-scale  bend¬ 
ing  of  the  (broken-up)  beam.  This  suggests  that  there  is  more  energy  in 

certain  (anisotropic)  turbulence  scales  than  any  of  the  spectral  models  pre¬ 
dict. 


*The  small-wavenumber  dependence  of  the  modified  von  Karmann  is 
(*  +  l/i>o2)  11/6,  while  that  of  Ref.  25  is  K~11/6^K  +  1/h  j  -11/6. 

these  expressions,  fc  =  L  /2tt . 


25.  D.  P.  Greenwood  and  D.  0.  Tarazano,"A  Proposed  Form  for  the  Atmospheric 
Microtemperature  Spatial  Spectrum  in  the  Input  Range ".RADC-TR-74-19 
February  1974.  * 
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There  is  clearly  a  need  for  further  analytical  work  and  experimentatioT 
on  the  problem  of  outer-scale  effects,  which  will  be  important  for  infrared 
systems  operating  over  reasonable  pathleng^hs. 

Inner  Scale  Effects 

Inner-scale  (i  )  effects  are  governed  in  the  Lutomirski-Yura 
3  4  u 

formulation  ’  by  an  expression  similar  to  that  for  the  outer  scale.  A  new 

critical  length  is  decined,  and  the  inner  scale  effects  are  determined 

by 


(40) 


which  is  identical  in  form  to  Eq.  (38).  In  order  to  avoid  the  effects,  we 

require  z/z^<  1,  which  at  visible  wavelengths  and  our  pathlength  of  1.6  km 

may  not  always  be  true.  Hence,  depending  upon  the  size  of  the  inner  scale, 

some  of  the  data  presented  below  may  have  been  affected.  However,  as  p  is 

reduced  due  to  stronger  and  more  highly-developed  turbulence,  Jt  will  tend 

o 

to  decrease  also. 

Retro-Reflector  Effects 

26 

It  has  been  reported  that  image  dancing  or  "angular  scintilla¬ 
tion"  from  a  laser-illuminated  retro-reflector  are  substantially  greater 
than  is  observed  from  a  nearby  beacon.  For  that  reason,  the  wander-tracking 
results  presented  below  were  obtained  with  a  point-beacon  which  is  coincident 
with  the  target  point-receiver. 

We  intend  to  perform  quantitative  experiments  on  the  retro  case, 
and  we  point  out  here  that  R.  Lutomirski27  has  preliminarily  analyzed  the 
effect  for  the  near-field  case.  It  appears  that  this  is  an  important  phe¬ 
nomenon  for  retro-tracking  systems. 


?6.  J.  P.  Hansen  and  S.  Madhu,  Applied  Optics  11,  233-238,  February  1972. 
27.  R.  Lutomirski,  private  communication. 
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B.  Experimental  Results 

In  this  section,  we  discuss  experimental  results  for  mean  irradiance 
and  fading,  with  and  without  tracking,  obtained  with  a  15  cm  transmitter 

over  a  1.6  km  path  at  6328  X.  The  range  of  the  independent  parameter  (D/p  ) 

<  <  o 

is  1  ^  (d/pq)  %  where  the  actual  variable  is  pq  or  strength  of  turbu¬ 

lence.  As  discussed  below,  results  for  smaller  values  of  this  parameter 
require  a  shorter  range,  which  has  been  established  and  is  currently  being 
utilized.  We  also  discuss  certain  aspects  of  the  wander  signal  per  se. 

As  a  matter  of  practical  necessity,  our  actual  transmitter  consists 
of  a  truncated-gaussian  with  an  aperture  diameter  of  15  cm.  The  experimental 
results  reported  here  are  in  preliminary  form,  in  that  we  have  not  yet 
determined  the  quantitative  effect  of  t  lis  truncation.  An  effective  diameter 
can  certainly  be  found  for  the  first  moment  of  irradiance;  a  different  diam¬ 
eter  will  apply  to  the  second  moment. 

For  the  sake  of  initial  simplicity,  the  diameter  D  utilized  in 
plotting  the  data  vs.  D/pQ  was  taken  as  the  total  aperture  diameter  of  15  cm. 
It  is  apparent  that  the  effective  diameter  will  be  less  than  this  value,  so 
that  such  data  points  will  need  to  be  moved  to  the  left.  This  will  be  seen 
to  improve  the  theoretical-experimental  comparison  in  some  cases. 

1  •  Real-Time  Irradiance  Signals 

The  smoothed  on-axis  irradiance  vs.  time,  with  and  without 
wander-tracking,  is  shown  for  the  weak  turbulence  case  ^  •  0.97^  in  Fig. 

10.  In  the  non-tracked  case,  th*  beam  is  initially  centered  in  terms  of  its 
wander  excursions,  and  It  gradually  drifts  off  due  primarily  to  slow  atmos¬ 
pheric  refractive  effects.  The  unsmoothed  irradiance  is  shown  in  Fig.  11a, 
and  the  corresponding  (ac-coupled)  logarithmic  signals  in  Fig.  lib.  The  low- 
frequency  nature  of  the  wander  component  of  the  signal  fluctuations  is  evi¬ 
dent. 

Similarly,  smoothed  signals  for  the  strong-turbulence  case 
<D/po  “  63)  are  shown  in  Fig.  12.  The  unsmoothed  linear  and  logarithmic 
signals  are  shown  in  Figs.  13a  and  13b  respectively.  It  is  apparent  that, 
in  the  high-turbulence  case  (beam  break-up),  the  advantage  of  tracking  is 
primarily  limited  to  the  longer-term  refractive  wander. 
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Figure  11a.  IJnsmoothed  irradiance  vs.  time,  for  the  case  of  Fig.  10  (lower  trace 
is  tracked).  The  abscissa  is  0.5  sec/div. 
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Fig. 


lib.  Log  irradiance  vs.  time,  for  the  ease  of  Fig.  10  (lower  trace  is 

Ordinate:  One  decade/div. 

Abscissa:  0.5  see/div. 


t  rac ked  ) . 
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(1-sec)  Ir rad  lance  vs.  time,  with  (lower  trace)  and  without 
racking  (upper  trace),  for  strong  turbulence. 


Figure  13a.  Unsmoothed  irradiance  vs.  time,  for  the  case  of  Fig.  12  (lower  trace  is 
tracked).  The  abscissa  is  0.5  sec/div. 
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Figure  13b. 


Log  irradiance  vs.  time 


Ord inate : 
Absc issa : 


for  the  case  of  Fig.  12  (lower  trace  is  tracked). 

one  decade/div. 

0.5  sec/div. 
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2.  Mean  I r radiance 


The  mean  on-axis  irradiance  for  the  non— tracked  and  tracked  cases 
is  shown  in  Fig.  14  as  a  function  of  (D/p^).  These  values  were  corrected 
for  transmitter  output-power  variations  from  run-to-run.  Since  the  aperture 
D  was  tixed,  the  appropriate  theoretical  curves  for  comparison  are  those  of 
^NT,T  norma^^ze^  by  Iq  (Figs.  1  and  3).  As  discussed  below,  small  values  of 
D/pq  could  not  be  achieved  over  this  path  with  reasonable  turbulence,  and 
hence  IQ  is  not  accurately  known.  However,  from  Fig.  3,  we  may  expect  that  the 
Largest  IT  at  small  (D/pq)  is  approximately  equal  to  T  and  this  normalization 
is  used  for  fitting  the  ordinate  to  the  data.  The  previously-mentioned  bias 
in  the  abscissa  is  evident.  It  is  noted  that  wander-tracking  has  more  of  an 
advantage  than  predicted  at  high  D/p^,  as  is  generally  observed  in  these 
experiments.  This  also  agrees  with  the  observations  of  Gilmartin  and  Holtz?8 
Also,  a  portion  of  the  1^  data  points  may  be  depressed  somewhat  owing  to 


longer-term  refractive  effects,  i.e.,  which  may  show  up  in  the  2  minute  aver¬ 
ages  used.  This  is  discussed  further  below. 

-2  — 

The  dependence  of  I  for  increasing  (D/p  )  is  generally  borne 


o  -  °  '-""o 

out.  For  lower  strengths  of  turbulence,  the  dependence  of  instantaneous 

beam  spread  ( -v  IT  *)  has  been  empirically  stated2^  to  be 
.  _  r  0.85  0.62,  0.65  CJ 

^  Y  =  2  k  ;  we  find  a  correlation  coefficient  of  0.89  for 


n 


y  <  10,  thus  corroborating  that  observation.  For  y  >  10,  the  correlation 

coefficient  is  6.8  x  10  ,  but  we  expect  an  approach  to  a  (p  dependence 

in  that  regime.  The  correlation  coefficient  for  ( p  _2  \  is  0.59  for 

\  NT  o  / 

—  >  10. 
p0 

In  Fig.  15,  we  show  the  ratio  ^  ^j^NT )  vs’  For  comparison,  the 

theoretical  curve  of  Fig.  4  is  also  shown.  The  tracking  advantage  is  gener¬ 
ally  higher  than  expected;  more  data  are  being  taken. 


3.  Fluctuations  in  Irradiance 

The  normalized  variance  of  irradiance,  with  and  without  tracking, 
is  shown  vs.  (D/ p Q )  in  Fig*  16a.  These  results  compare  nicely  with  the 
curves  of  Fig.  9,  which  are  shown  again.  For  comparison  with  the  non- 
tracked  case,  we  also  show  the  curve  from  Fig.  5  (Ref.  16). 


28.  T.  J.  Gilmartin  and  J.  Z.  Holtz, "Focused  Beam  and  Atmospheric  Coherence 
Measurements  at  10. 6p  and  0.63p  *  preprint,  Lincoln  Laboratories, 
October  1973. 

J.  A.  Dowling  and  P.  M.  Livingston,  J.  Opt.  Soc.  Am.  63,  846  (1973). 
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Figure  14.  Experimenta 


Figure  15.  Experimental  values 


D/p 

2  o 

Figure  16a.  Normalized  irradiance  variance  with  (o  ) (•)  and  without  (0/  )  (A)  tracking,  vs.  D/p  .  The  theoretical 
curves  from  Figs.  5  ( - )  and  9  T  ( - )  are  NT  also  shown.  0 


The  ratios  of  the  corresponding  tracked  and  non-tracked  variances  of 
Fig.  16a  are  shown  in  Fig.  16b.  The  advantage  remains  substantial  for 

—  >>1,  which  again  suggests  large-scale  wander  effects  greater  than  are 
included  in  the  theory.  For  the  tracked  case,  the  asymptote  is  seen  to  be 
unity  as  expected;  this  would  also  be  true  for  the  non-tracked  case  at  suf¬ 
ficiently  large  (D/p^).  More  data  are  being  taken. 

It  may  be  pointed  out  that  the  reduction  in  variance  for  the  tracking 
case  is  much  more  substantial  if  the  target  signal  is  low-pass  filtered,1  or 
alternatively,  if  a  finite  (but  small)  target-receiver  aperture  is  used. 

If  we  process  the  fading  signal  directly  for  the  log  amplitude,  rather 

than  the  linear  irradiance,  and  then  replot  the  abscissa  in  terms  of  the 

2 

Rytov  log-amplitude  variance  (o  )  for  a  point  source,  we  obtain  the  curve  of 

Xt  2  -3 

Fig.  17.  In  this  graph.  D/p  -  1  corresponds  to  o  -  3.7  x  10  ,  and  for 

°  XT 

values  substantially  larger  than  this  we  expect  log  normality  and  a  simple 

relationship  between  the  ordinates  of  Figs.  16a  and  17  (i.e.,  Eq.  31  with 

2  2 

B  “1).  In  particular,  the  asymptote  o^.  *  1  corresponds  to  o  *  0.17, 

X 

and  2  corresponds  to  0.27.  It  is  interesting  to  note  that,  for  the  non- 

tracked  case,  which  is  not  log-normal  at  small  values  of  the  abscissa  in 

Fig.  17,  the  value  of  the  log  amplitude  variance  is  more  or  less  constant 

regardless  of  the  turbulence  level.  In  such  a  case,  we  expect  the  same  to 

be  true  for  the  dynamic  range  of  the  fading. 

The  probability  distributions  for  log  irradiance  are  shown  in  Fig.  18 

for  the  tracking  and  non-tracking  cases  in  weak  turbulence  (D/pQ  =  0.97). 

The  residual  non-log-normality  in  the  tracking  case  is  related  to  the  details 

of  fading  for  an  approximately-diffraction-limited  focal  spot.  The  curve 

for  the  non-tracking  case  remains  to  be  interpreted  in  terms  of  the  statis- 

13 

tics  of  wander-fading.  For  strong  turbulence  (D/pq  *  173),  nearly-log- 
normal  distributions  are  observed,  especially  with  wander-tracking  (Fig.  19). 

The  spectra  of  the  log  fading  signal  for  weak  turbulence  (D/p  =  0.97) 
are  shown  in  Fig.  20,  for  the  tracking  and  non-tracking  cases.*  The  signifi¬ 
cant  spectral  components  are  below  e.g.,  10  Hz,  and  are  sharply  reduced  by 
tracking,  indicating  that  wander  is  predominant.  It  is  also  interesting  to 


*In  the  spectra  of  Figs.  20-22,  the  transverse  wind  speed  was  not  sufficient¬ 
ly  large  to  measure  meaningfully,  while  in  Fig.  23,  it  was  1.4  m/sec.  The 
general  conclusions  drawn  here  are  not  dependent  on  wind-speed  normalization. 
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Figure  19.  Cumulative  probability  (C.P.)  for  the  log  irradiance  in  strong  turbulence, 
for  the  wander-tracked  ( - )  and  non-t racked  ( - )  cases. 
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note  that  at  some  (fractional-Hz)  frequency,  there  is  an  apparent  peak  in  the 
non-tracked  spectrum,  coincident  with  a  dip  in  the  tracked  soectrum.  This 
feature  is  noted  in  the  other  examples  to  follow,  and  is  thus  related  to  a 
peak  in  the  spectrum  relating  to  wander  alone.  The  same  curves,  normalized 
by  the  log  amplitude  variance  (a  )  for  each  case,  are  shown  in  Fig.  20a. 

A 

For  this  low-turbulence  case,  the  spectral-crossover  between  wander-  and 
scintillation-predominance  is  evident  but  not  dramatic.  In  Fig.  20b,  we  show 
the  curves  of  Fig.  20  with  the  ordinate  weighted  by  the  frequency;  the  ten¬ 
dency  toward  separate  peaks  for  wander  and  scintillation  effects  is  apparent 
in  the  non-tracking  curve. 

In  Figs.  21-23,  we  show  the  log  fading  spectra  for  progressively 
increasing  turbulence  (D/pQ  =5.7,  20,  and  63  respectively).  As  beam  break¬ 
up  becomes  more  severe,  so  that  scintillations  predominate  over  wander,  the 
spectral  pairs  tend  to  flatten  out  and  come  together,  although  the  low-fre¬ 
quency  advantage  from  wander-cancellation  remains.  In  addition,  as  the  tur¬ 
bulence  grows,  there  is  a  tendency  for  the  spectra  of  the  tracked  and  non- 
tracked  cases  to  cross-over  at  higher  frequencies;  this  may  be  due  to  tracker 
response  to  beam-break-up  scintillations  and  lack  of  total  correspondence 
between  the  beacon  image  at  the  tracker  and  the  actual  distribution  of  energy 
at  the  target.  This  in  turn  corresponds  to  isoplanatism  effects  as  they 
relate  to  the  present  application  of  reciprocity,  and  will  be  discussed  fur¬ 
ther  in  a  future  report.  The  normalized  and  frequency-weighted  spectra  cor¬ 
responding  to  Fig.  23  are  given  in  Figs.  23a  and  23b  respectively  (rf.  Figs. 

20a, b).  With  strong  turbulence,  the  double  peak  is  no  longer  evident. 

4 .  Wander  Characteristics 

The  control  signal  in  the  wander-tracking  system  is  an  analog  of 
beam  wander.  An  example  of  a  real-time  horizontal  wander  signal  is  shown  in 
Fig.  24a.  In  some  cases,  the  wander  is  strongly  anisotropic,  which  can  be  a 
manifestation  of  both  slow  vertical  refractive  effects,  and  of  large-scale 
Lq)  turbulence  anisotropies.  An  extreme  example  of  the  former  is  shown 
in  Fig.  24b ,  which  is  from  the  same  data  run  as  Fig.  24a.  This  extreme  is 

typical  of  night-time  operation,  but  was  observed  to  a  much  lesser  extent  in 

some  of  the  daytime  runs  which  comprise  most  of  the  data  of  this  report. 

More  typically,  the  rms  vertical  wander  was  less  than  the  horizontal  value, 
which  is  probably  related  to  the  anisotropy  of  turbulence  at  larger  scales. 
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Figure  20a.  Curves  of  Fig.  20,  normalized  by  total  log  amplitude  variance  for  each 
case. 
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Figure  2 3n .  Curves  of  Fig.  23,  normalized  by  total  log  nmpli 
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Figure  24b.  Vertical  wander  signal  vs.  time  for  a  night-time  run  in  low  wind. 


The  slow,  refractive  components,  which  have  a  period  of  half  a  minute  or 
more,  are  difficult  to  characterize  in  terms  of  their  relationship  to  turbu¬ 
lence  effects  vs.  changing  thermal  stratifications.  The  effect  does  not  ap¬ 
pear  in  our  fading  data,  in  which  0.1  Hz  high-pass  filtering  was  effectively 
used.  However,  our  measurements  of  I  were  taken  by  centering  the  open- 
loop  transmitter  on  the  target  axis  and  then  averaging  for  120  sec.;  in  this 
case,  slow  wander  can  reduce  the  apparent  mean  irradiance  over  that  which 
would  be  obtained  with  shorter  averaging  times,  and  therefore  enhance  the 
ratio  Qualitative  examination  of  real-time  irradiance  records  does 

not  indicate  that  this  effect  was  large,  for  120  sec  averaging  times,  and 
shorter  intervals  results  oc  course  in  greater  data-spread.  Ideally,  each 
measurement  of  1^  would  be  repeat  d  a  number  of  times,  using  shorter  aver¬ 
aging  times  and  re-centering  ea._h  Lime. 

A  quantitative  distinction  between  the  two  vertical  beam-bending 
mechanisms  can  be  made  using  the  probability  distributions  of  the  wander 
signal.  In  Fig.  25a,  we  show  the  distribution  for  the  horizontal  signal  of 
Fig.  24a,  with  equivalent  high-pass  filtering  at  four  different  cut-off  fre¬ 
quencies.  The  distributions  are  gausslan,  with  coefficients  of  skewness  and 

-1  -4 

kurtosis  for  the  100  sec  case  of  3.6  x  10  and  3.1,  respectively.  In 
Fig.  25b,  we  show  the  distributions  for  the  vertical  signal  of  Fig.  24b, 
where  the  non-gaussian  nature  is  quite  apparent.  The  best  approximation  to 
a  gausslan  distribution  is  seen  to  occur  for  an  inverse  cut-off  frequency 
of  20  seconds. 

A  summary  of  the  characteristics  of  the  wander  signals,  averaged  over 

a  number  of  runs,  is  shown  in  Table  I.  The  ratio  of  average  horizontal/ 

vertical  (rms)  wander  angles  Is  1.6.  The  correlation  coefficient  refers  to 

the  basic  theoretical  prediction  (wander  term  of  Eq.  20),  namely  (mean-square 
2 

wander  angle)  ^  .  In  spite  of  the  Influence  of  slow  refractive  effects, 

the  vertical  wander  exhibits  much  less  bias  than  does  the  horizontal.  This 
may  be  due  to  variable-wind  effects,  which  influence  the  horizontal  case 
only. 

The  average  two-dimensional  mean-sqjare  wander  angle,  for  the  runs 
included  in  Table  I,  implies  a  value  of  | |  In  Eq.  (20)  of  3.69,  which  is 
nearly  equal  to  that  deduced  in  the  semiquantitative,  phenomenological 
reasoning  of  a  previous  section. 
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TABLE  1.  AVERAGE  CHARACTERISTICS  OF  BEAM  WANDER  SIGNALS 


Correlation  Coefficient 

Skewness  Kurtosis  (Exper.  vs.  Theory) 

Vertical  2.0  x  10~2  2.76  0.90 

Horizontal  2.3  x  10~2  2.89  0.36 

Total  0.97 


The  power  spectrum  and  frequency-weighted  spectrum  of  the  wander  signal 
of  Fig.  24a are  shown  in  Fig.  26.  The  theoretical  phase-difference  spectrum, 
which  is  closely  related,  is  predicted  to  fall  off  with  an  exponent  of 
(-2/3, -8/3)  for  frequencies  above  breakpoints  (f  )  respectively: 


f  n,  \ 
1  L 

o 


f„  *  \ 

D 


(41) 


where  v^  is  the  perpendicular  component  of  wind  velocity.  For  frequencies 

below  f  ,  the  behavior  is  very  sensitive  to  the  detailed  behavior  of  the 

wind.  Since  this  was  a  low-wind  case,  we  expect  a  log-log  slope  of  -8/3; 

the  actual  slope  was  -2.78  with  a  correlation  coefficient  of  0.98.  In  Fig. 

27,  horizontal  and  vertical  spectra  are  shown  for  a  higher  wind  (f0  ** 

^horizontal 

4  Hz),  and  the  -2/3  slope  is  well  supported  for  the  horizontal  case.  Tne 
vertical  'ase  shows  a  greater  slope,  and  we  note  that  f£  for  the  vertical 
wind  is  approximately  zero  near  the  ground.  In  other  cases,  such  as  shown 
in  Fig.  28  ( f 2  ^  0) »  both  the  horizontal  and  vertical  speccra  fall  between 
the  two  theoretical  slopes. 

5.  Transmitter  Focus  Effects 

As  discussed  in  a  preceding  section,  the  first-order  theory  pre¬ 
dicts  that  a  focused,  near-field  transmitter  with  wander-cancellation  will 
result  in  smoothing  of  scintillations.  It  is  now  understood  taat  this 
effect,  if  it  occurs,  can  only  be  realized  for  the  condition  D/p  <  1.  In 

_ _ _  O 

30.  A.  J.  Huber,  "Measurements  of  the  Temporal  Power  Spectra  of  a  Propagated 
10.6  micron  Wavefront",  Technical  Report,  Rome  Air  Development  Center, 
Spring  1973. 
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Figure  26.  Power  spectrum  S(f)  ( - 

for  wander  signal  of  Fig 
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Figure  27.  Power  spectrum  of  horizontal  ( - )  and  vertical  ( - )  wander  signals 

for  significant  horizontal  wind.  ( . )  indicates  f “^^dependency . 


-63- 


F  iguri  28. 
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for  low  wind  speed,  exhibiting  slopes  intermediate  to  the 
values  of  -2/3  ( . )  and  -8/3 
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fact,  we  expect  that  the  actual  requirement  is  —  <<  1,  so  that  phase  dis- 

Po 

tortion  is  negligible.  As  a  follow-up  to  our  worn  as  reported  in  Ref.  22, 

where  qualitative  evidence  of  such  smoothing  was  seen,  we  are  conducting 

experiments  on  thf  (wander-tracked)  signal  fading  vs.  transmitter  focus 

adjustment.  As  discussed  in  that  reference,  the  smoothing  effect  is  expected 

to  be  very  critically  dependent  upon  focus  adjustment. 

Attempts  to  observe  this  effect  are  shown  in  Figs.  29a  and  b.  The 

parameter  a^z  is  simply  z/R,  where  R  is  the  radius  of  curvature  of  the 

transmitter  wavefront.  Since  (D/pQ)  was  on  the  order  of  3.5  for  these  runs, 

the  effect  was  inconclusively  observed. 

As  discussed  below,  we  are  now  operating  on  a  much  shorter  path.  As 

seen  in  Fig.  30  a,b,  the  predicted  effect  is  observed,  including  the  large 

increase  in  scintillations  with  slight  mis-ad justments.  However,  we  have  not 

yet  operated  with  —  «1,  and  we  have  not  observed  a  reduction  to  a  value  of 
2  po 

o  substantially  below  that  for  a  point  source  on  the  same  path. 

XP 

The  point-source  scintillations  are  necessarily  quite  low  under  these 
conditions,  and  any  residual  fading  due  to,  e.g.,  tracker  error  or  jitter 
will  be  quite  important.  How  *r,  the  tracker  resolution  is  0.02  x  the  dif¬ 
fraction  limit,  and  we  still  hope  to  see  substantial  smoothing.  This  effect 

1  7_1  Q 

has  important  practical  implications  for  uplink  beams,  especially  at 

longer  wavelengths. 


6.  Further  Experiments  and  Parameter  Variations 

More  data  are  obviously  required,  especially  for  (D  £  p  ).  How¬ 
ever,  the  interesting  case  is  that  of  large  Fresnel  number,  D>>  (z/k)1^2. 

We  thus  require  that  the  following  ratio  be  large: 


✓  z/k 


0.411 


2 

Point-Source 

Rytov 


1.44  k 


_7 

10 


11 


n 


3 

5»1  .  (42) 


2 

In  order  to  achieve  this  at  a  reasonable  Cn  (i.e.,  reasonably  well- developed 

turbulence),  we  require  long  wavelengths  and  short  pathlengths,  with  the 

stronger  dependence  on  the  latter.  If  we  specify  p  =  D  =  15  cm,  we  can 

2  ° 

calculate  the  corresponding  C  vs.  psthlength  and  wavelength,  and  the 

l/2n  1/2 

resultant  value  of  (2irz/k)  =  (A„  .  This  is  shown  in  Table  II. 
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Experimental  log  amplitude  variances  (o  )  divided  by  the  first-order  theoretical  values  for  a  point 
source  (a  2),  Vs.  focus  condition  anz  *E  (Ref.  22).  The  value  of  D/p  is  3.5. 


TABLE  2. 

REQUIRED  C  VS.  RANGE  AND 
n  WITH  RESULTANT 

WAVELENGTH  FOR  COHERENCE 
FRESNEL  NUMBER 

SCALE  =  15  CM 

z 

A 

c2 

n 

(A,)1/2 

1.6  km 

6328  A 

2.7  x  10'j* 

9.9  x  io':r 

7.7  x  10'U 

3.2  cm 

1.6  km 

3.8m 

7.8  cm 

1.6  km 

10.6m 

13  cm 

100  m 

6328  A 

4. A  x  lO'Jij 

i.6  x  io':, 

1.2  x  10 

8.0  mm 

100  m 

•  3.8m 

1.9  cm 

100  m 

10.6m 

3.2  cm 

It  is  clear  that,  in  order  to  achieve  good  (pq > D)  data  with  a  15  cm 
transmitter  and  suitable  Fresnel  number,  the  combination  of  a  longer  wave¬ 
length  and  short  pathlength  is  desirable.  We  are  currently  utilizing  a  91  m 
path  at  6328  A,  and  we  will  subsequently  modify  the  tracking  transmitter  to 
accommodate  10.6m.  These  pathlength  and  wavelength  variations  will  also 

serve  to  extend  our  comparisons  with  theory,  which  have  previously  been 

2 

limited  to  variations  in  C  . 

n 

In  accordance  with  Eq.  (AO),  since  the  cases  of  interest  now  involve 
Po'  cm»  inncr  scale  effects  should  be  negligible.  At  10. 6p  for  the 
short  path,  pq  may  become  substantially  larger  than  15  cm,  permitting  data 
for  the  pQ  »  D  case;  however,  according  to  Eq.  (38),  outer  scale  effects 
may  then  become  noticeable.  This  problem  may  be  reduced  by  decreasing  D 
and  accepting  a  lover  Fresnel  number. 

C.  Future  Work 

In  the  above  we  have  described  future  experimental  work  with  a 
shorter  pathlength  and  longer  wavelength,  and  further  theoretical  efforts 
on  the  second  moment  of  irradiance  as  well  as  the  details  of  wander-tracking 
with  a  gaussian  beam.  In  addition,  we  intend  to  determine  the  scintillation 
characteristics  from  generalized,  extended  sources  such  as  laser-illuminated 
diffuse  reflectors.  This  will  include  the  strength  and  (transverse)  struc¬ 
ture  of  scintillations,  which  are  important  for  remote-probing  applications 
with  non-cooperative  targets. 

As  a  related  item,  we  plan  to  return  to  the  problem  of  saturation 
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1  31 

of  scintillations  over  a  long  path,  *  using  very  small  receivers  and 

large  (e.g.,  10  kHz)  bandwidths.  We  will  compare  the  results,  including 

32  33 

covariances,  with  recent  physically-based  analyses  by  Yura  and  Clifford, 

and  we  expect  to  snow  that  a  new  similarity  theory  appearing  in  the  Russian 
34 

literature  is  incorrect.  The  results  are  again  important  for  remote¬ 
probing  applications. 

Finally,  we  will  return  to  such  problems  as  outer  scale  effects,  angular 
scintillation  of  retro  reflectors,  and  isoplanatism  as  it  relates  to  point- 
target  reciprocity  and  wander-tracking. 


31.  J.  R.  Dunphy  and  J.  R.  Kerr,  J.  Opt.  Soc.  Am.  62,  981-986,  August  1973. 

32.  H.  T.  Yura,  J.  Opt.  Soc.  Am.  64_,  59-67,  January  1974. 

33.  S.  F.  Clifford,  G.  R.  Ochs,  and  R.  S.  Lawrence,  J.  Opt.  Soc.  Am.  64, 
148-]54,  February  1974. 

34.  M.  E.  Gracheva,  et  al ," Similarity  Correlations  and  Their  Experimental 
Verification  in  the  Case  of  Strong  Intensity  Fluctuations  of  the  Laser 
Radiation  ".preprint,  August  1973. 
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HI*  SHORT-TERM  MICROi^ERMAL  AND  SCINTILLATION  STATISTICS 

Since  the  preceding  report  on  this  program,1  considerable  progress  has 
been  made  on  the  problem  of  short-term  microthermal  and  scintillation  statis¬ 
tics  and  their  interrelationship.  The  analytical  considerations  and  experi¬ 
mental  data  will  appear  in  the  next  report,  including  data  to  be  taken  during 
the  forthcoming  summer.  In  this  section,  -e  will  briefly  review  the  progress 
which  has  been  made. 

Confidence  Intervals 

We  can  now  analytically  relate  the  data-spread  or  confidence  intervals 
which  are  obtained  from  short-time-average  measurements  of  strengths  of  tur¬ 
bulence  and  scintillation  respectively.  In  particular,  this  involves  the 
characterization  of  the  spatial  field  of  mean-square  microthermal  fluctua¬ 
tions.  Significart  experimental  data  have  been  obtained. 

Averaging  Time  Effects 

We  have  resolved  the  issue  of  the  apparent  inconsistencies  in  ave  aging- 
time  theory  as  applied  to  measurements  on  intermittent  turbulence  and  result¬ 
ing  scintillations.  In  addition,  we  can  quantitatively  specify  the  degree 
of  "intermittency".  Fourth  moments  of  probability  distributions  are  neces¬ 
sarily  involved,  as  are  the  integral  scales  of  squared  and/or  smoothed  ran¬ 
dom  processes. 

Computer  Simulation 

We  have  formulated  a  means  of  computer  simulating  ensemble  members  of 
the  instantaneous  propagation  path,  within  the  short-path  or  geometrical- 
optics  realm.  These  simulations  will  be  performed,  and  related  to  the  con¬ 
fidence-interval  considerations  mentioned  above.  In  addition,  attempts  will 
be  made  to  include  the  time-variation  (over  a  fixed  period)  as  part  of  the 
definition  of  each  ensemble  member,  and  finally,  to  extend  the  technique  to 
the  first-order  physical-optics  realm. 
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IV.  PUBLICATIONS  AND  PRESENTATIONS 


Forthcoming  publications  and  presentations  on  work  related  to  this  pro 
gram  are  as  follows: 

1.  James  R.  Dunphy  and  T.  R.  Kerr,  "Atmospheric  Beam  Wander  Cancella 
tion  by  a  Fast  Tracking  Transmitter",  to  appear  in  Journal  of  the 
Optical  Society  of  America,  June,  1974. 

2.  J.  R.  Kerr,  "Mean  Irradiance  and  Fading  for  Near-Field  Laser 
Transmitters  with  Beam  Wander  Cancellation:  Unified  Analytical 
Treatment",  Pape-  TuA2,  Topical  Meeting  on  Optical  Propagation 
Through  Turbulence  (sponsored  by  the  Optical  Society  of  America), 
July  9-11,  1974,  Boulder,  Colorado. 

3  J.  R.  Dunphy  and  J.  R.  Kerr,  "Mean  Irradiance  and  Fading  for  a 
Focused,  Near-Field  Laser  Transmitter  with  Beam  Wander  Cancella¬ 
tion:  Experimental  Results",  Paper  ThA3,  ebove  conference. 

4.  Philip  Pincus,  R.  A.  Elliott,  and  J.  R.  Kerr,  "Short-Term  Sta¬ 
tistics  of  Turbulence  and  Optical  Propagation",  Paper  WB8,  above 
conference. 

In  addition  to  the  above,  written  papers  on  each  of  these  areas  are 
contemplated  in  the  fairly  near  future. 
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